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^. . Abstract 

We present results from numerical simulations of the (2 + l)-dimensional Gross-Neveu model 
with a f/(l) chiral symmetry and Nf = 4 fermion species at non-zero temperature. We 
provide evidence that there are two different chirally symmetric phases, one critical and one 
with finite correlation length, separated by a Berezinskii-Kosterlitz-Thouless transition. We 
have also identified a regime above the critical temperature in which the fermions acquire a 
screening mass even in the absence of chiral symmetry breaking, analogous to the pseudogap 
behaviour observed in cuprate superconductors. 



1 Introduction 

Phase transitions in QCD and model field theories have been studied intensively over the 
last decade both analytically and numerically. Understanding the phase diagram of QCD is 
becoming increasingly important in view of recent experimental efforts to create and detect 
the quark-gluon plasma in the relativistic heavy ion collisions at BNL and CERN. 

Since the problem of chiral symmetry breaking and its restoration is intrinsically non- 
perturbative the number of available techniques is limited and most of our knowledge about 
the phenomenon comes from lattice simulations. Because of the complexity of QCD with 
dynamical fermions, studies so far have been done on lattices with modest size and in various 
cases the results are distorted by finite size and discretization effects. The three-dimensional 
Gross-Neveu model (GNM3) has been proved to be an interesting and tractable model to 
study chiral phase transitions both numerically by means of lattice simulations and analyti- 
cally in the form of large- iVj expansions [1—10]. 

In this letter we present results of numerical simulations of the ^(A^j)v-invariant GNM3 
with a f/(l) chiral symmetry at non-zero temperature. This model is described by the 
following continuum Lagrangian density (we work in Euclidean space throughout): 

£ = ^,(^ + mo)vI/, - ^[(^,,vI/,)2 _ (^,75iI/,)2]. (1) 

We treat \E'j, ^j as four-component Dirac spinors and the index i runs over Nj fermion 
species. In the chiral limit mo— >0, the U{1) invariance is given by 

^^e^^^^^ ; ^^^e'^^\ (2) 

This Lagrangian is considerably easier to treat, both numerically and analytically, if it is 
bosonized by introducing auxiliary fields a and tt: 



C = ^i(^ + mo + a + i75vr)^» + TT^l^^ + ^^)- (3) 



At tree level, the fields a and tt have no dynamics; they are truly auxiliary fields. However, 
they acquire dynamical content by dint of quantum effects arising from integrating out the 
fermions. The model is renormalizable in the 1/Nf expansion unlike in the loop expansion 

i- 

Apart from the obvious numerical advantages of working with a relatively simple model 

in a reduced dimensionality there are several other motivations for studying such a model: (i) 
at T = for sufficiently strong coupling it exhibits spontaneous chiral symmetry breaking and 
the pion field vr is the associated Goldstone boson; (ii) the spectrum of excitations contains 
both baryons and mesons, i.e. the elementary fermions and the composite fermion— anti- 
fermion states; (iii) the model has an interacting continuum limit for a critical value of 
the coupling l/g"^, which has a numerical value ~ 1.0/a in the laige-Nj limit if a lattice 
regularisation is employed [^ ; (iv) numerical simulations of the model with chemical potential 
/i 7^ show qualitatively correct behaviour, unlike QCD simulations 0, ^. 

The phase diagram of GNM3 with various global symmetries at non-zero temperature 
and density has been studied extensively in [3—10]. More specifically [|], ^, it was shown 



that the thermally induced phase transition of the Z2-syninietric model belongs to the two- 
dimensional Ising universality class in accordance with the dimensional reduction scenario 
||Tl| which predicts that the long-range behaviour at the chiral phase transition is that of 
the {d — 1) spin model with the same symmetry, because the IR region of the system is 
dominated by the zero Matsubara mode of the bosonic field. Outside the non-trivial region 
the critical properties of the system are described by mean field theory according to the large- 
Nf prediction. It was also shown that the two-dimensional Ising scaling region is suppressed 
by a factor Nj'^^^ g. 

At leading order in 1/A^/ the U{1) model has a second order chiral phase transition at Tc = 
^^ |l[], where ruf is the fermion dynamical mass at zero temperature. The leading order 
effective potential has the same form as the discrete symmetry case with the replacement 
cr^ — s> a'^ + ir^. This conclusion is expected to be valid only when Nf is strictly infinite, 
i.e. when the fluctuations of the bosonic fields are neglected, since otherwise it runs foul of 



the Coleman- Mermin- Wagner (CMW) theorem |T^, which states that in two-dimensional 
systems the continuous chiral symmetry must be manifest for all T > 0. Next-to-leading order 
calculations demonstrated that the symmetry is restored for arbitrarily small T<0 [7—10]. 



As argued in [^ |T0| the model is expected in accordance with the dimensional reduction 



scenario to undergo a Berezinskii-Kosterlitz-Thouless (BKT) transition W3\ at Tbkt which 



is associated with the unbinding of vortices like in the two-dimensional XY model. In terms 
of the reduced temperature t = T — Tbkt the scaling behaviour of the correlation length, 
susceptibility and the specific heat is given |13[ by 



m-e'''-'\ x{t)--e-^ a ~r/^ + constant, (4) 

where for t -^ O"*", v = 1/2, t] = 1/4 and a = —dv = —1. It is easier to visualize this 
scenario if we use the "modulus-phase" parametrization a + in = pe*^. In two spatial 
dimensions logarithmically divergent infrared fluctuations do not allow the phase 6 to take 
a fixed direction and therefore prevent spontaneous symmetry breaking via {6) ^ 0. The 
critical temperature Tbkt is expected to separate two different chirally symmetric phases: a 
low T phase, which is characterized by power law phase correlations (^e^^i^'ie'^^^^^) ~ x'^^"^^ at 
distances x ^ 1/T but no long range order (i.e. a spinwave phase where chiral symmetry is 
"almost but not quite broken"), and a high T phase which is characterized by exponentially 
decaying phase correlations with no long range order. In other words for < T < Tbkt there 
is a line of critical points characterised by a continuously varying < rjiT) < i. However, 
since the amplitude p is neutral under f/(l), the dynamics of the low temperature phase 
do not preclude the generation of a fermion mass rrif oc p whose value may be comparable 
with the naive prediction of the large- A^/ approach; this was demonstrated in d = (1-1-1) 
in [Q, in which results from an exactly soluble fermionic model ||T^ are generalised to 
GNM2. The effect may be understood in terms of the spectral function p{s). For orthodox 
mass generation, the spectral function in the broken chiral symmetry phase is a simple pole 
p{s) oc 6{s — mV)- By contrast in the BKT scenario |T^ the function is modified to a branch 
cut p(s) oc {^s — mf)'^~^6{y/s — mf); although chiral symmetry is manifest, the propagating 
fermion constantly emits and absorbs massless scalars and hence has indefinite chirality. 



In [1^ it is shown that in GNM3 with T > 0, next-to-leading order corrections cause 



fermion mass generation to occur for T <T^, where T* > Tbkt, ie. it predicts a "pseudogap" 



phase in which the system is non-critical, chiral symmetry is manifest, and yet the energy 
required to create a spin-i excitation does not vanish. Such behaviour is interesting because 
it is characteristic of the underdoped phase of certain cuprate superconductors [1^ . 



We have attempted to check the vahdity of (a) the BKT scenario, and (b) mass generation 
without chiral symmetry breaking, by performing lattice simulations. Our results with Nf = 
4 show that the symmetry is restored for a large range of non-zero temperatures. The 
susceptibility of the order parameter diverges at a temperature which we associate with 



Tbkt- We also verify the prediction of [|I0[ that the fermions have non-zero mass in a region 
where the system is not critical, i.e. for T > Tbkt- Our simulations have not shown evidence 
for a second phase transition to a massless phase at T = T*, as predicted by large- A^/ methods 
10| , although we cannot rule it out. 



2 Numerical Simulations 

The model in its bosonised form can be formulated on the lattice using the following action: 

Slat = E xix)M,,yx{y) + - E x{x)x{x)[Y, a{x) + ie{x) ^ 7r(£)] + ^ E[^'(^) + ^'(^)]' 
^,y X (x,x) (x,x) 3 X 

(5) 

where x ^^^ X ci-re Grassmann-valued staggered fermion fields defined on the lattice sites, 

the auxiliary field a is defined on the dual lattice sites and the symbol (x, x) denotes the set 

of 8 dual lattice sites x surrounding the direct lattice site x . The fermion kinetic operator 

M is given by 



Mx,y = - 



^y,x+b "y,x-0 



+ nJ2 Vi^ix) [Sy,x+0 - Sy,x-o] , (6) 

^ !/ = l,2 



where r]^{x) are the Kawamoto-Smit phases (^—i^^o+-+x,^-i g^j^d the symbol e(x) denotes 
the alternating phase (— l)^o+a;i+a::2^ ^j^g simulations were performed by using the standard 
hybrid Monte Carlo algorithm in which complex bosonic pseudofermion fields $ are updated 
using the action ^^{M^M)~^^. According to the discussion in |^, simulation of A^ staggered 
fermions describes Nf = AN continuum species; the full symmetry of the lattice model in 
the continuum limit, however, is U{Nf/2)v ® U{Nf/2)v ® U{1) rather than U{Nf)v ® U{1). 
At non-zero lattice spacing the symmetry group is smaller: U{Nf/4)v (S> U{Nf/4)v ® U{1). 
In this study we have used N = 1. 

In order to study the behaviour of the chiral symmetry first at T = and then for 
T > we set the fermion bare mass to zero. Without the benefit of this symmetry breaking 
interaction, the direction of symmetry breaking changes over the course of the run so that 
E = yj^x^^ix) cind n = yJ2x'^{x) average to zero over the ensemble. It is in this way 
that the absence of spontaneous breaking of this continuous symmetry on a finite lattice is 
enforced. The next best thing to measure is an effective "order parameter" |$| = vSM-lP, 
which is a projection onto the direction of $° = (S, 11) separately for each configuration 
||T7| . In the case of chiral symmetry breaking |$| differs from the true order parameter (S) 
extrapolated to the chiral limit, because in the absence of a symmetry breaking term it 



is impossible to disentagle the fluctuations of the order parameter field from those of the 
Goldstone modes. Including the effects of the latter will give |$| > Sq, where Sq denotes 
the value of S in the chiral limit. The simulations at T = were performed on 24^ x 36 
lattices for (3 = l/g"^ = 0.45 - 0.775. In Fig.[^ we plot |$| and the fermion mass m/ versus 
p. We fitted these two observables to the scaling functions |$| = ai(/5^""^ — /5)^'""s and 
ruf = a2{(3^^^ - I3f for /5 = 0.625 - 0.750. We extract p = 1.08(5), /3„„g = 1.1(1) which 
are consistent with the values of the large- A^/ prediction /5 = z/ = 1 [Q. The values of 
the associated bulk critical couplings are: (3^'^''^ = 0.87(3) from the order parameter fit and 
pbuik _ Q_8g(^x) from the fermion mass fit. At very strong coupling (small f3) the data deviate 
slightly from the fitted scaling functions due to discretization effects. 

In order to study the effects of non-zero temperature we performed simulations on asym- 
metric lattices with constant temporal extent Lt = A and spatial extents Lg = 30, 50, 100, 150. 
We vary T by changing (3. At /3^""^ the lattice spacing becomes zero and T —>■ oo. Simulations 
of the Z2 model with Lt = 4 show a thermally induced phase transition at /5f^ = 0.565(3), 
which should be compared with the value jSc = 0.76 expected in the Nf —^ 00 limit on a 
lattice with Lt = 4. For Lg = 100 and 150 we accumulated approximately 30, 000 to 45, 000 
trajectories for the range of couplings (3 < 0.54 and 10, 000 to 20, 000 for larger values of 
j3. The trajectory length was approximately 1.5 and was chosen at random from a Poisson 
distribution in order to decrease autocorrelation times. For Lg = 30 and 150 we chose a 
trajectory length of approximately 1.0 and accumulated 40, 000 to 50, 000 trajectories for 
f3 < 0.54 and 20, 000 to 30, 000 for larger values of /3. In Fig.0 we plot |$| versus (3 for the 
various lattice sizes together with the results from simulations of the Z2 model on lattice 
sizes 4 X 50^ and 4 x 100^. It is clear that the order parameter of the Z2 model is independent 
of the lattice size until just before the transition at /5 = /5f^, whereas in the f/(l) model |$| 
has a strong size dependence for a large range of values of (3, i.e. it decreases rapidly as the 
spatial volume increases in accordance with the expectation that chiral symmetry should be 
restored for T > 0. We believe that the signal for chiral symmetry breaking at very low T 
is a result of finite size corrections. The finite spatial extent Lg provides a cut-off for the 
divergent correlation length and according to the BKT scenario the slow decay of the cor- 
relation function ^e«^(^)e-«^(o)^ with exponent i]{T) < 0.25 for T < Tbkt ensures a non-zero 
magnetization even in a system with very large size. In other words on small lattices the 
spinwave phase is expected to look like a broken phase. A similar effect has been observed 



in the two-dimensional Gross-Neveu model with a discrete symmetry at non-zero T [18|: For 
large- A^/ the density of kinks (whose condensation is responsible for the restoration of chiral 
symmetry at T > 0) is exponentially suppressed and therefore on lattices with spatial extent 
smaller than the kink size the system looks as if it is in the broken phase. 
In Fig.^ we plot the susceptibility of the order parameter 

measured on lattices with spatial size Lg = 30, 50, 100, 150. It is clear that x has a peak 
for Lg = 30,50,100 at PsKTiLg) ~ 0.54,0.53,0.515 respectively. The values of x iisar the 
Lg = 150 peak are very noisy and we don't plot them in this graph. This is clear evidence of 
a phase transition, with a critical coupling Pbkt significantly less than f3^'^. The transition 
occurs at a much lower temperature than the critical temperature of the GNM3 with a Z2 



symmetry, because the infrared fluctuations are stronger in the continuous symmetry case. 
Another interesting observation is that in the low temperature phase the susceptibility has 
stronger size dependence than in the high temperature phase and the error bars at low T 
are much larger that at high T despite the fact that the statistics at low T are much larger 
than at high T. This is evidence that the system is critical in the low T phase in accordance 
with the BKT scenario. Another useful quantity to measure is the so-called specific heat C^, 
which we calculated from the fluctuations of the bosonic action Sb = ^J2xW^{^) + 7r^(x)]. 
C^ is given by 

C. = ^{{S'b) - {S^n (8) 

and is plotted versus coupling for Lg = 50, 100, 150 in Fig.^. It is clear that (a) C^ has a 
broad peak at /3 ^ 0.50 and (b) it does not show any divergent behaviour or significant 
finite size effects which is consistent with the BKT scenario, according to which a = —2. 



Analogous behaviour has been observed in numerical studies of the XY model [|T9 . 

Now we discuss the issue of fermion mass generation. Since the temporal direction of our 
lattices is short, in order to study the asymptotic behaviour of the fermion propagator we 
have chosen to measure the spatial fermion correlator Gs{x) along one of the spatial lattice 
axes, which yields information on spatial correlations in the thermal medium rather than 
the spectrum of excitations. The quantity analogous to mass is the so-called screening mass 
Ms] for free, massless continuum fermions this quantity is given by the lowest Matsubara 
mode for fermions, ie. Mg = Mq°"* = ttT. On a finite temporal lattice this becomes Mq"* = 
sinh~^ [sin(7r/Lt)] . For Lt = A the screening mass in the absence of dynamical mass generation 
is Mg*^*^ 0.658. If a dynamical mass M is generated, then naively we expect 

Ml = Ml + M^. (9) 

We calculated the screening masses from the exponential decay of the spatial correlation 
functions using an ansatz for the fit which is motivated by the form of the free propagators 

Gs{x) = A{1 - (-1)") sinh(M,(x - LJ2)) + B{1 + (-1)") cosh(M,(x - LJ2)). (10) 

The screening masses extracted from the simulations of GNM3 with U{1) symmetry on 
lattices with varying Ls, together with the results of the Z2 model from simulations on a 
4 X 100^ lattice, are shown in Fig.^ It is clear that the fermions have non-zero screening mass 
M in the sense that Mg > Mq, even when the order parameter is very small. The screening 
mass is independent of Ls for almost all values of /3 - in particular we observe that M^'-^^ is 
comparable with M^^ for (3 < (3^"^. In the Z2 case we ascribe mass generation to orthodox 
chiral symmetry breaking. For U{1) the situation is more subtle, as shown in Fig.^ where 
we plot the spatial fermion correlator of both models on a 4x 100^ lattice at /3 = 0.51. The 
Z2 data are non-vanishing on all timeslices, consistent with the standard scenario of mass 
generation via broken chiral symmetry. By contrast, the vanishing of Gs{x)^^^'^ for even x, 
corresponding to i? = in (p!OD, signals a manifest chiral symmetry, which follows because 
the t/(l)e symmetry of staggered fermions implies that the only non-vanishing elements of 
the propagator are Geo and Goe- By comparing Fig.| with Fig.|^ we infer that the fermions 

6 



are massive, ie. M > 0, in the phase where both chiral symmetry is restored and the system 
is not critical, ie. the susceptibihty is finite and size independent. This imphes the existence 
of a pseudogap phase for jSbkt ~ /3 ~ Pc^ ■ 

Strictly speaking, in demonstrating that Mg > Mq we have not excluded the possibility 
that the relation (|) could be replaced by M^ ~ A{T)M^ + M'^ with A- 1 > M'^/M^, imply- 
ing a temperature-dependent modification of the speed of light rather than mass generation 
2T[. The similarity of M^^^^ and M^^ at low T where A ^ 1 suggests this is unlikely to 



be the dominant effect for (3 ^ P^^; however for /3 ^ 0.56 the two models begin to diverge, 
and the level ordering switches. Since the departure from the free-field result in this regime 
is greater for the U{1) model (which has twice as many scalar fields) as for the Z2 model, 
it is conceivable that the effect is due to in-medium modifications of the massless fermion 
propagator, and could be calculable in a large- A^/ framework. To understand the physics of 
this regime the two effects need to be disentangled via a more detailed study of the dispersion 
relation for the spatial correlator, ie. the value of Mg for several distinct Matsubara modes. 
This will require simulations closer to the continuum limit, i.e. on lattices with Lt > 8. 
For this reason we cannot at this stage identify a second phase transition at /5 = /5* where 
mass generation switches off. It is possible that such a transition may be an artifact of the 
large- A^y approach [n]|. 



3 Summary and Outlook 

Our numerical study of the f/(l)-symmetric GNM3 provided evidence for: (a) existence of 
two chirally symmetric phases separated by a transition which is possibly a BKT transition in 
agreement with the dimensional reduction scenario; (b) dynamical mass generation without 
symmetry breaking in accordance with the prediction of the large- iVj calculation and (c) the 
existence of a pseudogap phase where the system is not critical (i.e. T > Tbkt) and the 
fermions acquire a non-zero mass. Our results are another example of what was shown by 



Witten ||lj] in GNM2 with a U{1) symmetry at T = 0, i.e. that the large- iVj expansion is a 



reliable guide to the properties of the model as long as one interprets the results carefully. 



In [Q it was shown that although large- iV/ predicts spontaneous symmetry breaking which 
is in contradiction with the CMW theorem, the scalar field correlator falls off like x~^^^f , 
i.e. there is an "almost long range order" together with fermion mass generation which does 
not break chiral symmetry because the physical fermion is a superposition of positive and 
negative chirality states and has zero net chirality. The results of |^, where it was shown 
that for GNM3 with a Z2 chiral symmetry at non-zero T the region (in the space of m/T) 

— 1/2 

of the large- A^/ behaviour squeezes out the two-dimensional Ising region by a factor Nr , 
also support the spirit of Witten's statement. 

We are currently extending our work in various directions. We wish 

• to increase the statistics near and below the transition in order to perform finite size 
scaling to extract the exponents ri{T). 

• to repeat the simulations with L^ = 8 in order to study the fermion dispersion relation. 



to study the A^/-dependence of our various results, which can be compared with analytic 
predictions |[T( 



We also plan to study the GNM3 with an SU{2)l ® SU{2)ii chiral symmetry at non-zero 
temperature. According to the dimensional reduction scenario this model's critical properties 
at non-zero T should be the same as the two-dimensional spin model with an 0(4) symmetry, 
i.e. it is expected to have a transition at T = 0^, and to have dynamical mass generation 
for T > 0. 
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Figure 1: Order parameter |<l>| and fermion mass rrif at T = vs. coupling j3 on a 24^ x 36 
lattice. 
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Figure 2: Order parameter |$| vs. f3 for L^ = 4 for both U{1) and Z2 symmetric models. 
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Figure 3: Susceptibility x vs. p. 
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Figure 4: Specific heat Cy vs. /?. 
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Figure 5: Fermion screening mass Mg vs. (3. The horizontal hne shows the lowest Matsubara 
mode M^"*. 
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Figure 6: Fermion spatial correlators Gs{x) extracted from simulations of U{1)- and Z2- 
symmetric GNM3 on 4 x 100^ lattices at /? = 0.51. 
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